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Abstract 

We provide a general formalism to calculate the infrared correlators of multiple 
interacting scalar fields in the de Sitter space by means of the stochastic approach. 
These scalar fields are treated as test fields and hence our result is applicable to the 
models such as the curvaton scenario where the fields that yield initially isocurvature 
modes do not contribute to the cosmic energy density during inflationary expansion. 
The stochastic formalism combined with the argument of conformal invariance of the 
correlators reflecting the de Sitter isometries allows us to fix the form and amplitude 
of the three-point functions completely and partially for the four-point functions in 
terms of calculable quantities. It turns out that naive scaling argument employed 
in the previous literature does not necessarily hold and we derive the necessary and 
sufficient condition for the correlator to obey the naive scaling. We also find that 
correlation functions can in principle exhibit more complicated structure than argued 
in the literature. 



1 Introduction 



De Sitter spacetime is as fundamental as Minkowski spacetime. It describes the accel- 
erated expansion of the Universe sourced by the cosmological constant, or the effective 
vacuum energy mimicked by some scalar field. In addition to the observation that present 
Universe is undergoing accelerated expansion, it has become a standard paradigm that the 
small inhomogeneity, i.e., tiny deviation from the Friedmann-Lemaitre- Robertson- Walker 
Universe, is also thought to be explained by conversion of the quantum fluctuations of 
light scalar fields generated during de Sitter expansion in the early Universe PQ. Because 
of this, it is important to study statistical properties of the scalar field in de Sitter space. 

The simplest scenario for the generation of the curvature perturbation is to assume 
that inflaton which causes inflation simultaneously generates observed amplitude of the 
primordial curvature perturbation ~ 10~ 5 pp. Despite its simplicity and that its predic- 
tions are consistent with observations j2], this scenario is not a prediction of (unknown) 
fundamental theory but rather an assumption. Actually, there are also many alternatives 
to this scenario most of which introduce other scalar fields which are originally isocurva- 
ture modes and convert to the curvature perturbation after inflation. The curvaton [3H5] 
and the modulated reheating models [6j[7] are the representative ones that belong to this 
category. The scalar fields in this category are generically negligible during inflation, that 
is, they do not affect the dynamics of the inflationary expansion. Because of this, the 
correlation functions among such scalar fields become invariant under the de Sitter isome- 
tries. Note that this is not necessarily true for the case of the inflaton fluctuation since 
the inflaton can affect the expansion of the Universe. 

The purpose of this paper is to study generic statistical properties of the correlation 
functions of scalar fields which enjoy full de Sitter invariance. In particular, we are con- 
cerned with infrared limit of the correlators. This is equivalent to large distance limit or to 
late time limit in de Sitter space since any two points that are close originally is eventually 
stretched to arbitrary large distance by accelerated expansion. 

Some literature [SHTTj already discuss the generic shape of the correlators consistent 
with de Sitter invariance. In these papers, the scalar fields appearing in the correlators 
are (implicitly) assumed to scale like (—r]) Aa for —rj <C H^ 1 (late time), where t] is the 
conformal time of the de Sitter metric given by 

and A a is a constant specific to the scalar field <fi a and cannot be constrained only by the 
de Sitter invariance. For instance, (0 a ( ? 7i, a ? i)0b( ? ?2, ^2)) ~ ( — Vi) Aa (~ 7 72) Afc etc., where A a 
neither depends on the number of fields nor the type of fields appearing in the correlator 
under consideration but just depends on </>„. This time dependence of the correlators then 
allows us to (partially) fix their dependence on the spatial coordinates by requiring that 
the correlators are invariant under the de Sitter isometric transformations. However, it is 
not clear how wide class of models satisfy the above scaling. 
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In order to evaluate the correlation functions for the de Sitter invariant state without 
introducing the scaling assumption a priori, we make use of the stochastic formalism. This 
formalism was introduced and developed in [T2HT6] . It has been since then employed for 
various inflationary models to study the infrared behavior of the scalar fields (e.g., [T&H25] ). 
The formalism has been also used to determine the distribution of the initial value of 
the scalar field that may become important in the late Universe. This is important, for 
example, to provide the initial value of the quintessence field which causes the accelerated 
expansion of the present Universe [26|[27]. 

The stochastic formalism solves the infrared dynamics by treating the long wavelength 
modes as the classical statistical variables that are sourced by stochastic noises coming 
from the short wavelength quantum modes. Since the dynamics can be solved, we can fix 
the form of the correlation functions unlike in the literature where they are constrained 
only by the argument of the de Sitter invariance. The stochastic formalism is especially 
useful when the nonperturbative effect is crucial to obtain the correct correlation functions 
as in the case of massless self-interacting scalar field [TB]. This is simply because infrared 
dynamics is solved without invoking perturbative expansion, that is, nonperturbative effect 
is automatically taken into account in the stochastic formalism. 

Despite the many existing applications of the stochastic formalism to particular in- 
flationary models, we do not find any paper that discusses general consequences of the 
stochastic approach on the infrared properties of the correlators of the multiple interact- 
ing scalar fields, which motivated us to address this issue. This article is the report of 
the calculations of the two-, three- and four-point functions of the multiple scalar fields 
derived by the use of the stochastic formalism. As will be demonstrated explicitly, in 
addition to the scaling index of the correlators, their amplitudes, which are completely 
unconstrained within the framework of the symmetry argument, can also be expressed in 
terms of quantities that are reasonably calculable in the stochastic approach. Thus we can 
completely fix the (infrared) correlation functions. We will find that the naive universal 
scaling <p a ~ (— r/) Aa in the correlators does not always hold and provide the necessary 
condition for this scaling to hold by explicitly constructing three-point and four-point 
functions. 

2 Basics of the stochastic formalism 

Our purpose is to study superhorizon evolution of (weakly) interacting multiple fields in de 
Sitter space and to calculate the resulting correlation functions in the stochastic formalism. 
Before doing this, in this section, let us briefly review the basic points of this formalism. 
For more details, see, for example, [T2T[T6] . 

The first step of the stochastic formalism is to split the field (operator) into the long 
wavelength part <j) L and the short wavelength part; 



<j>(t,x) = <j> L (t,x) + 



1 
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T 6(k - ea(t)H) [a^(t)e 
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+ h.c. 
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where 6{x) is a step function and we have used a bold letter = (0i, </> 2 , • • • , 0tv) to make 
it refer to the multiple scalar fields. The second term on the right-hand side represents the 
contribution from short wavelength modes whose wavenumbers are greater than ea(t)H. 
Since we are assuming that the background spacetime is de Sitter one (a(t) = e ), any 
mode that originally belongs to the short wavelength part eventually enters <j) L . Although 
the evolution of the short wavelength part depends on the nature of the interactions, in 
this paper, we assume that any interaction is negligible for the short wavelength part and 
<j)£ obeys the massless Klein-Gordon equation in de Sitter space whose solution is given by 



a,k 



y/2k\' k 



where the subscript a of runs from 1 to N and rj = — is the conformal time. Then 
we can treat and at as the standard annihilation and creation operators. 

The second step is to consider the evolution equation for the long wavelength part 
<p L . Applying the slow-roll approximation and neglecting the higher-order terms in the 
short-wave modes and mode mixing terms, it reads 

Mt,x) = -^V a ( ( / ) ) + f a (t,x), (4) 

where V a = dV/d<fr a , and we have abbreviated the subscript L for the coarse-grained field. 
This equation may be regarded as a classical Langevin equation with a stochastic noise 
term, f a , which is given by time derivative of the short wavelength part and represents 
modes of k = ea(t)H that join <j) L at time t. Straightforward calculation shows that f a is 
random Gaussian whose two-point function is given by 

H 2 

(f a {ti,xx)f b {t 2 ,X2)) = Sat—ditt - t 2 )j (ea(t)H\xi - x 2 \), (5) 

where jo(x) = is the unnormalized sine function. 

Then we can easily show that the one-point probability density pi((j)(x),t) = p(0i(x), ^(x), - ■ ■ ,t) 
obeys the Focker-Planck equation: 

Introducing the dimensionless potential by v(<f>) = 47r 2 V(<f>)/ (3H 4 ), a general solution of 
the above equation can be written as 

pi(4>, t) = e~ v ^ «n$n(0)e- A " (i - io) , (7) 

n 

where $ n (0) is an eigenfunction of A r -dimensional Schrodinger equation: 

( 1 d 2 I, ,\ ... 4vr 2 A n ^ ... 
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To be definite, we only consider the case the eigenfunctions are normalizable and satisfy 



d(j) $ m (»$„(» = 5 mn . (9) 
The left-hand side of Eq. ([8]) can be written as 

?K-^ + "»)(^ +t '-)*" w - (io) 

By multiplying $„(0) to the above expression from left and integrating it by parts, the 



left-hand side of (IE]) becomes the integral of 



2 



. Therefore, A„ > 0. In 



particular, the eigenfunction $o having the minimum eigenvalue (i.e., Aq = 0) is given by 

%((f>) =Me- v ^\ M = (J ' dxj) e ~ 2v ^ \ (11) 

At sufficiently late time, all the modes having the positive A n decays in Eq. ((7|) and p\ 
becomes independent of time; 

Pl (0,t)->p cq (0)=ArV 2 ^. (12) 

Here we used = M which holds from the normalization J d/j)pi((j>, t) = 1. This is a 
distribution function for an equilibrium state achieved in the de Sitter space. 

It is possible that the integral appearing in the definition of diverges. In such a 
case, Eq. ([7]) does not possess any static solution and hence there is no equilibrium state. 
This happens, for example, for massless free scalar field. 



3 Correlation functions 

The formalism explained above enables us to evaluate correlation functions between mul- 
tiple fields. In the following, we will derive the infrared behaviors of the two-point, three- 
point and four-point functions separately by means of the stochastic formalism. Derivation 
of the expression of the two-point functions mostly follows the one developed in [16], which 
also showed their de Sitter invariance constructed from the equilibrium distribution func- 
tion. 

3.1 Two- point functions 

The (de Sitter invariant) spatial correlators (on superhorizon scales) at equal time can be 
written as 

{<l>a(t,xi)<l> b (t,x 2 )) = / d4 l d4> 2 4*14)1 p 2 (</>\</> 2 ,£), (13) 
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where p2(0 1 ,$ 2 ,t) is the probability density of finding (j) 1 at Xj (i — 1,2). It can be shown 
that p 2 obeys the following equation, 



i=l 



-P3 



# 3 9 2 P3 



:j (ea(t)^|xi -f 2 |). (14) 



At early times when the points x\ and $2 are deeply inside the same Hubble patch, j on 
the right-hand side becomes unity. For such a case, it can be shown that 

p 2 (0 1 ,#,t)=5(0 2 -0 1 ) Pcq (0 1 ), (15) 

constitutes a static solution. This is very reasonable since setting jo = 1 means that <j>i 
and <f>2 are fully correlated, i.e., <j>i = <j> 2 and the appearance of p eq reflects Eq. ffl2l that 
guarantees the de Sitter invariance [H]. Therefore, Eq. (|15p can be used as the initial 
condition for Eq. 014|) . 

It is hard to find the analytic solution of Eq. (Tl4|) with the initial condition (fl5l) . Here 
we make an approximation that 

jo(ea(t)Hr) ~ 9(1 - ea{t)Hr). (16) 

This drastically simplifies the equation without losing any essential point of the stochastic 
formalism. This approximation allows us to write down the solution of Eq. f[T41) after x,\ 
and X2 are separated by super-horizon distance; 

p 2 ((/) 1 ,<f> 2 ,t) = J cZ0 r n(^ 1 ,t;^ r .,t r )n(0 2 ,t;0 r ,t r .)p eq (^ r ), (17) 

where t r is a solution of ea(t r )H \xi — x 2 | = 1 and represents the time when at X\ and (j> 
at X2 get uncorrelated. Here n^ 1 , t\\ <f> 2 , t 2 ) is the transition probability from <j) = <j> 2 at 
t = t2 to <j> = (f) 1 at t = t\. Its expression in terms of the eigenfunctions is given by [T6] 



00 

2^ 



H(0U;#,t 2 ) = e-^H"^) ^$ n (0 1 )$ n (0 2 )e" A "^ 



-ta) 



n=0 



This expression of p 2 appeals to our intuition, that is, p 2 is given by the product of the 
probability of <f> going to <f>\ from <j) r by the stochastic process described by Eq. (j3j) and 
that of going to , from r with a weight p eq (0 r ). Using this picture for p 2 , we have 



(0a(*,£l)&(*,£ 2 )) = y #'# 2 0a0' y d^M^^t-^^t^U^^t-^^t^p^). (19) 

Substituting Eq. (TIB]) into Eq. (HPJ, we find 

(Mt,x x )Mt,^)) ^Arf2A^At l) (HR 12 )- 2 ^ L ex V , (20) 

n=0 ^ ' 

5 



where R12 = a(t)\xi — x 2 \ is the physical distance between X\ and x 2 and A a is defined 
by 

A™ = j d<j> 4> a e- v ^ n (<f>). (21) 

Equation (120]) is our expression for the two-point functions. To minimize the effect of e, we 
choose e so that it satisfies exp (— ^-lne) ~ 1 for the dominant mode contributing to the 
two-point functions, as suggested in [16]. Now let us consider correlation function ( 12 Oj) on 
sufficiently late time or (equivalently) large scales in which case a(t)\xi —X2\/H~ 1 is quite 
large. Then the leading contribution to the correlator is from a state $n(</>) labeled by 
an integer n having minimum A n (apart from the ground state n = which has A n = 0) 
with nonvanishing Aa A b ■ Generally, n depends on the choice of fields (0 a , and can 
vary for different set of fields of the correlators. In particular, it may happen that n for 
A a A^ which we denote by n a is different from m and A n , a 7^ A nb . In such a case, n for 
A a n ^A^ which we denote by n a b may or may not exist. If n a b does not exist, it means 
no correlation between <j) a and 0&. Alternatively, our result can also allow a possibility 
that n a f, exists, in which case <j) a and <pb are correlated. By a simple consideration, we find 
that such n ao is either equal to or larger than n a or n& whichever is greater. In short, two 
different fields having different n (i.e., n a ^ rib) can in principle have correlation between 
them, which is consistent with our findings for the simple example demonstrated in the 
Introduction. 

On the contrary, if some of the fields are interacting so that the minimum integers for 
the nonvanishing of Aa for such fields are all the same, then the correlator exhibits a 
universal behavior in the sense that the scaling index for any field takes the same value 
and is completely given by A^, eigenvalue of the A r -dimensional Schrodinger equation. In 
this case, correlators at late time scales as 

(Mt,xi)Mt,X2)) ~ Wi-x 2 \ H ■ (22) 



3.2 Three-point functions 

What we want to evaluate is the spatial three-point functions evaluated at equal time t; 

(Mt,xi)Mt,^)<i>c{t,x 3 )) = J # x # 2 # 3 <f> x Ml p 3 (0 1 ,0 2 > 3 ,O, (23) 

where psfy 1 ,^ ,^ 3 ,t) is the probability density of finding <j>' at Xi (i = 1,2,3). It can be 
shown that p% obeys the following equation, 

-ET=y. -EjT ~^TT PZ ) o 2 P>MP>M ~A 2 ait a jj J °( 6a W H \ X i- X 3\)- 1 24 ) 

dt ^ \ 31 J 8tt 2 d(p a d<j)l\ 4vr 2 j^j d^dcpl 

At early times, all the three points are inside the Hubble radius, i.e., a(t)\xj — Xj\H 1 
for any 1 < i, j < 3, and the fields are maximally correlated each other. During this epoch, 
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the evolution equation for p% can therefore be well approximated by Eq. (|24|) with all the 
jo being replaced by unity. We can verify that this equation allows the following solution, 

p,{<j>\<j> 2 )( j>\t) = 5(<t> 3 -<t> 2 )5(<f> 2 -<f> 1 )p C(l (<t> 1 ), (25) 

whose physical meaning is obvious from the reasoning we made earlier. This solution is 
independent of time and can be used as a de Sitter invariant initial condition of Eq. (1241) . 
Then the problem is to solve Eq. (1241) with such an initial condition until sufficiently late 
time when all the points are separated by super-horizon length and any correlation between 
the fields at different points is turned off. Although this is a well defined mathematical 
problem and we can in principle solve Eq. (12^1) and perform the integrals appearing in 
Eq. (123]) to get the three-point functions, we find it difficult in practice to solve Eq. (T21|) 
which is highly involved partial differential equation. Fortunately, as long as we are only 
concerned with sufficiently late time behavior of the three-point functions, which is ac- 
tually the present case, there is a way to derive an analytic expression without directly 
solving Eq. f|2^|) as we will demonstrate below. The point is to utilize the de Sitter isome- 
tries which allows us to find the three-point functions at general points by implementing 
the coordinate transformation that preserves the de Sitter metric from some extreme con- 
figuration of the points where analytic evaluation of the three-point functions (to a very 
good approximation) is feasible. 

As is well known, there are 10 isometries for the metric of the de Sitter spacetime. 
Among 10 isometries, translations and rotations for the spatial coordinate constitute 6 
isometries. We also have the dilatation isometry which amounts to multiply both time 
and spatial coordinates by the same constant factor. The remaining 3 isometries are 
complex mixing whose infinitesimal form is given by 

T]' — Tj — 2r](b ■ x), x' = x + b(—r] 2 + x 2 ) — 2{b ■ x)x, (26) 

where x 2 = x 2 and b is an infinitesimal constant vector. On super-horizon scales, or on 
sufficiently late time, in which r\ is much smaller than x, the finite version of the above 
transformation can be written as 

rf = Jl , f = * + X% . (27) 

1 + 26 ■ x + b 2 x 2 ' 1 + 2b ■ x + b 2 x 2 

The transformation of the spatial coordinate does not involve time and becomes exactly 
what is known as the special conformal transformation. The special conformal transfor- 
mation combined with the dilatation, rotation and translation transformations constitute 
the conformal transformation [28]. We will come back to this point later when we utilize 
the conformal symmetry to fix the correlators. 

Now let us consider the three-point function in the squeezed limit with different time 
coordinates for the different points; 

(MtuyiWbfamWbfam))- (28) 
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Although what we are interested in is the equal time correlators (t\ = t 2 = ^3 — t), for the 
moment, we let them to be independent due to the reason which will become clear later. 
To be definite, we take 2/3 = and |j7i| ^> I2/2I (squeezed limit), which is always possible 
without a loss of generality. For convenience, we define R = \y\\ and r = \y 2 \- All the time 
coordinates t\,t 2 and £3 are assumed to be very large so that any two different points are 
eventually separated by super-horizon size distance. Instead of directly solving Eq. (I24p . 
the following physical consideration enables us to evaluate Eq. (128]) . By definition, y 2 
and 2/3 (this is actually 0) are close together compared to y-±. Therefore, a field at 2/1 first 
gets uncorrelated and starts to evolve independently when the physical distance between 
y\ and other two points becomes equal to (eH)' 1 . Strictly speaking, the epoch when 
aH\yi — y 2 \ = 1 occurs is different from that when aH\yi — y 3 \ = 1 is satisfied. But the 
difference between these little affects the final result in the squeezed limit and we can 
safely take them as being equal. At this time, y 2 and 2/3 are still deeply inside the Hubble 
radius and fields at those two points take the same value. As the Universe expands, the 
physical distance between y 2 and 2/3 then becomes equal to (eH) -1 . After this time, all the 
fields at different points evolve separately. This picture, which is a good approximation 
when there is a huge hierarchy among the lengths of the sides of the triangle, enables us 
to write the three-point function in the following form, 



(Mti,yi)Uh,m)Mh,m)) = j # 1 # 2 # 3 # R # r <\>\<g<\>\ xi{<t>\t 1] <l> R ,t R )xi{<t>\t 2] <l> r ,t r ) 

xn(0 3 , t 3 ; <i> r , t r )n(4> r , t r - <j> R , t R ) Pcq (<t> R ). (29) 

Substituting Eq. (fl8~j) . the above expression reduces to 

£,m,n 

Xe -A<(*i-*H) e -A m (t 2 -t r ) e -A„(t s -t r ) e -A e (t r ~t R ) /^Q\ 

where Be mn is defined by 

B tmn = j ' d4> e v M<t> e (<j>)$ m (<j>)® n (<{>), (31) 

and is totally symmetric under the permutation of the indices. 

Now let us consider the late time behavior of Eq. f )30|) . Since all the eigenvalues satisfy 
A n > 0, the leading contribution comes from terms with a particular set of (£, fh, n) having 

nonvanishing Aa 'Bg^ with the minimum decay rate e~ A ^ 1 e~ A ™* 2 e~ A "* 3 . This 

means that each integer of (£, rh,n) is determined by the lowest value of the eigenvalues 
with the condition that A^ A^ A^ Bi mn does not vanish. Note that this condition does 
not necessarily fix (I, fh, n) uniquely apart from the trivial permutations and, depending 
on the interactions among the scalar fields, it is possible that there are more than one set 
of (£,fh,h). Generally speaking, all the numbers can be different from each other, can be 
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partially equal or completely coincide and concrete values of (£, in, n) needs specification 
of the underlying model. This may be understood by considering the simplest case, i.e., 
single field case in which Bm and Bn 2 vanish and the lowest contributions are either 
B 122 A^ A (2) A (2) + perms, (when A 3 > 2A 2 - Ai) or B 113 A W A w A i3) + perms, (when 
A3 < 2A2 — Ai). If the model yields A3 = 2A 2 — Ai by chance, both two contributions 
decay in time at the same rate and none of the two terms can be neglected even at 
sufficiently late time. We will come back to the single field case later. 
Using the equations for t r and £r given by 

t r = -^-\n(eRH), t R = --!- In (erH) , (32) 
ti ri 

we find that the three-point functions for sufficiently late time become 

(Mh, yi)Mt 2 , m)Mh, m)) « A/" 2 £ B Mn A^A^e A ^ +A ^ +A ^ 

min 

x(HR)- — (Hr) (33) 

where the summation indicated by "min" is done for all the possible sets of (£, m, n) 
satisfying the condition mentioned above. Here, we fixed e in the same way as the case of 
two-point functions. This is the late time three-point functions in the squeezed limit. 

This form allows us to write the expressions for the general shape of the triangle 
formed by x±, x 2 and x 3 by using that the left-hand side of Eq. (1331) is invariant under 
the transformation that preserves de Sitter isometry. To understand this, notice that 
we can always move one of the point, say x±, to a vector having very long length by 
performing the transformation (|27l) . Indeed, if we choose b as — % + £, then we find that 

the transformed point is y\ = £/£ 2 , whose distance from the origin can be arbitrary large 
in the limit £ <C 1. By this transformation, time coordinates, which have the same value 
in the original frame, take different values in the new frame. Thus, equal time correlator 
for arbitrary configuration of points is related to the squeezed correlator with different 
time coordinates for different points by isometry-preserving transformation. This is the 
reason why only the information in the squeezed limit is enough to obtain the three-point 
functions for any configuration of points. 

The de Sitter isometries for the spatial coordinates at late time \rj\ — > become the 
conformal transformation which serves a base for the conformal field theory [9]. In the 
language of the conformal field theory [2B], focusing on any one particular term in Eq. (1331 . 
the field cj) a can be interpreted as a primary field of a conformal weight —A a /H. Since it 
is well established how to obtain the general expression of the three-point functions of the 
conformal fields out of the squeezed limit [28], we do not expand the detailed discussion 
here and we only give the final result for the three-point functions which is given by 

(Mt,2i)Mt,&)4>c(t&)) « M 2 Y,BMnAfAf ) Af\HR 12 )- !l£l ^ 



x(HR 23 ) h ( H R 31 ) * , (34) 
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where Rij = a(t)\x*i — Xj\ represents the physical distance between Xi and Xj. This equation 
is one of our primary result. This shows that all the information regarding the three- 
point functions can be obtained once we know the eigenfunctions and eigenvalues of the 
N dimensional Schrodinger equation (jSJ). As a consistency check, we can verify that 
squeezed limit of Eq. (Ru = R13 = R, R23 = r) gives back Eq. ( 1331) . Also, we can 
implement the similar derivation of the three-point functions for the equilateral case in 
which \yi— 2/2 1 — W2— 1 = 1 2/3 — yi\- This case also allows the evaluation of the three-point 
function without resorting to the direct computation of Eq. (T24"l) . It can be verified that 
the result coincides with the equilateral case of Eq. ( |34|) . 

Now there are several points to be remarked. First, as mentioned earlier, (£,m,n) is 
given by the condition that it is a set of integers as small as possible with nonvanishing 
Aa^ A^ Bfrhn- I n principle, this integer set can vary for different choice of fields (a, b, c). 
This suggests that only a knowledge of correlators of the three product of the same field 
(i.e., (0 a 0a0a) etc.) is not enough to know the scaling behavior of the three-point function 
of the three different fields. Furthermore, if at least one of the field appearing in the 
correlator is different from others, it can happen that the correlator cannot be given by a 
single term with a power-law form. Instead the correlator becomes a sum of up to six terms 
each of which exhibits the different power-law behavior. Notice that all of those terms 
are not necessarily nonvanishing and it is also possible that only some of them remain 
nonzero. 

Secondly, the result (134j) is obtained without the use of the perturbative expansion 
in terms of the strength of the interactions among fields. In the standard approach, the 
three-point functions (and higher order functions as well) are calculated by using the so- 
called in-in formalism which usually uses perturbative expansion and truncation of the 
calculations at some order to yield an analytic expression. Although this approach is com- 
pletely justified as long as the higher-order terms contribute much less to the final result, 
it is known that some particular model (for example, massless scalar field with a quartic 
self-interaction) requires nonperturbative treatment to obtain the reliable correlators. In 
more general terms, the perturbative approach fails when the system does possess the de 
Sitter invariant state only if the interactions are present. In such a case, the standard 
perturbative approach needs some care, if not impossible, to get the correct result. On the 
other hand, our result (I34p is derived by the stochastic approach. As is well known, the 
stochastic approach includes the nonperturbative effects coming from the long wavelength 
modes. Except for some simple models, numerical computations are required to solve 
Eq. (JSJ) in order to obtain the eigenvalues A n . But whichever computation method is used, 
the obtained A n contains the nonperturbative effect. This is also true for the amplitude 
of the correlator, i.e., Aa A^ A^ B^n- The point is that everything is reasonably cal- 
culable in the stochastic approach while it is hard to take into account the higher-order 
or nonperturbative effects in the standard perturbative expansion of the in-in formalism. 
Therefore, our result will be quite useful when the nonperturbative effect is crucial to get 
the correct correlators. 
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Figure 1: Double squeezed quadrangle for which R u ~ R\ 2 ^> R23 — R24 ^ -R34, where 
Rij denotes the physical distance between $ and yj. 

3.3 Four-point functions 

Contrary to the case of the three-point function where consideration in the squeezed limit 
is enough to get the correlator for the general triangle, this is not the case for the four-point 
function. The reason behind this is that the transformation specified by the vector b is not 
sufficient to convert any quadrangle to the one having hierarchy among all the sides of the 
quadrangle. Thus we need to solve the evolution equation for the probability distribution 
function for four variables even to obtain the late time behaviour of the correlator. Yet, 
it would be interesting to see to what extent we can restrict the form of the four-point 
function in analytic way. 

Let us first consider the double squeezed quadrangle for which i? 14 ~ i? 12 3> -R23 3> -R34, 
where R^ denotes the distance between and yj (see Fig. [[]). For this quadrangle, using 
the similar reasoning we made in the case of the three-point function, the four-point 
function can be written as 



where X = i? 14 , Y = i? 2 4 and Z = R u (X ^> Y ^> Z). Interestingly, the four-point 
function in this limit is written solely in terms of the quantities characterizing the three- 
point functions. At sufficiently late time, this reduces to 




e -Ap(t 3 -tz)-Aq{t4-t z )-A a (t z -t Y )-A m (tY-tx) 



(35) 




xe 



— Apt3— Aqt4 



(HX) 




H 



(HZ) 



H 



(36) 
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where the meaning of the "min" in the summation is the same as the case for the three- 
point function. Now, the expression of the right-hand side of the above equation must 
be the double squeezed limit of the correlator for the general quadrangle. Using again 
the procedure of restricting the form of the four-point function in the conformal field 
theory [28], we find that the four-point function ( 1361) for the general quadrangle becomes 

(rW-)) =^E&(^.^)n^)^. (37) 

\ i=l I min x ' i<j 

where A x = ^p, A 2 = ^ etc. and A = A x + A 2 + A 3 + A 4 . The function /™"ff 304 
depending on the two arguments that are invariant under the de Sitter isometric transfor- 
mation (only for late time or on superhorizon scales) is completely arbitrary function at 
the level of the symmetry argument. Some more nontrivial information is needed to (even 
partially) fix the form of it. 

Although we cannot find analytic form of /™™ff 3a4 for the whole domain of the argu- 
ments, we can derive its asymptotic behavior or a value at specific point for some limiting 
cases of the arguments by considering the corresponding squeezed shape of the quadrangle. 
For instance, comparison between Eqs. ( 1361) and f|37|) leads to the following asymptotic 
form of f™%g 3ai {z, z) when z < 1 as 

■CXkuC*. *) * B^B^A^A^AjSA^-^^, for z « 1. (38) 

In a similar way, consideration of the other quadrangles with different shapes shown in 
Fig. [2] leads to the following expression; 

ffhnpq I -t\ ^ ,(m)4(n)4(p) 4(9) -^- + Am+A ^+' 

J 01020304^' / — J -'mns J ->spq^ ai - n - a2 ^a 3 ^o 4 * 

rihnpq M \ ^ P_ __P___ ^(n)4(p) 4(g) -^ + Am+A ^ 

J 01020304 \ > 41 1 — ±J mqs- L -'snq^-a- L ^o 2 ^a 3 ^a 4 ^ 

pfhnpq -i\ ^ rp A (m) A (n) A (p) A (q) 

J 01020304 I A > L J — - L mnpq- n -a 1 ^a 2 ^a 3 ^a 4 > 

In Eq. (HTj) . we have introduced T mnP q defined by 

T m npq = J d(j> e 2v ^ $ m (<f>) $ n (<j>) % (<j>) % (<f>) , (42) 

which is a new quantity to characterize the four-point function. These expressions for 
the several limiting cases indicate that four-point function having the double squeezed 
quadrangle is determined by B mns , i.e., related to (square of) the three-point function 
while the correlator having the single squeezed one is determined by T mnpq which has 
nothing to do with the three-point function. Finally, let us remark that Eqs. ( 138]) - ( 14"T|) 
cover all the possible limiting cases which we can obtain. One can show that we cannot 
make such a quadrangle which provides /™^ 304 (z, w) with z — 1, w <C 1, or z -C 1, w — 1 
or z, w 1. 



for ^ > 1, (39) 

for ^ > 1, (40) 
(41) 
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Figure 2: Left quadrangle having i? 12 ~ -R34 > R13 ~ -R24 yields Eq. (1391) . Middle 
quadrangle having Ru ~ -R34 ^> Ru ~ -R23 yields Eq. ( HUl) . Right quadrangle having 
i? 12 « i? 13 « R u > i? 23 = i?3 4 = R 24 yields Eq. gTJ). 

4 Discussion and summary 

In this paper, we have developed a general formalism to calculate three- and four-point 
functions among test light scalar fields in de Sitter space in the framework of the stochastic 
formalism. Stochastic approach treats the coarse grained fields as classical variables that 
are affected by the random Gaussian noises reflecting the effect that the short wavelength 
modes enter the long wavelength modes due to the accelerated expansion. We investi- 
gated behaviors of those correlators at sufficiently late time when the fields at different 
spatial points evolve independently and derived general formulae of the correlators. In the 
stochastic approach, one generally needs to solve the Focker-Planck equation to evaluate 
the time evolution of the correlator since its solution (probability density) appears in the 
definition of the correlator. This program (with a reasonable approximation) works for 
the two-point functions and the result is given by Eq. ( 1201) . This result clearly shows 
that the amplitude of the two-point functions is determined by defined by Eq. (12 lj) 
which represents the expectation value of 4> a for the n-th state e~ v ^Q n ((f>) corresponding 
to the eigenvalue A n which determines the scaling index. It is important to notice that by 
construction the stochastic formalism incorporates the dynamics of the long wavelength 
modes without expanding the basic equations in terms of the amplitude of the coarse 
grained fields or the strength of the interactions. Therefore, once Aa and A n and hence 
the two-point functions are obtained, they are the results that do not rely on the pertur- 
bative expansion and automatically involve the nonperturbative effects. Except for some 
simple cases, Eq. (JSJ) does not allow the analytic solution and numerical computation is 
generally required to obtain the concrete values of A^ and A n . However, solving Eq. (jSJ) 
is a mathematically well-defined problem and, in principle, it can be solved (especially 
for the case with a small number of fields) without any principal difficulty This plausible 
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feature of the stochastic formalism is in striking contrast to the standard calculation of the 
correlation function that usually uses perturbative expansion in terms of the strength of 
the interactions and truncates the expansion at some order. In most cases, the truncation 
is done at the lowest order that yields nonvanishing contribution and extension to the 
higher order entails very messy expression that requires careful consideration to obtain 
physically meaningful results. 

Direct manipulation of the evolution equation for the probability density becomes very 
difficult for the case of the three-point functions. Therefore, in this paper, we adopted 
another approach which works as long as we are only interested in the late time behavior of 
the correlators. The basic idea is to utilize the de Sitter isometries and to consider only the 
squeezed case in which one of the three spatial points is separated far way than the others. 
This limiting case allows us to build up the analytic expressions of the correlators at late 
time. We then converted the correlators to the ones for the arbitrary configuration of points 
by using the fact that the value of the correlator itself remains the same under the isometric 
transformation for the de Sitter invariant state. Among the ten de Sitter isometries, the 
three reduces to the special conformal transformation for the spatial transformation. It is 
this transformation that enables us to transform any triangle into the arbitrarily squeezed 
one and to obtain the analytic form of the (only late time) three-point function for any 
configuration of points which is given by Eq. ( l34l) . 

Now it would be interesting to consider the consequences of Eq. ( 1341) focusing on the 
single field case for simplicity. Even in this case, we find nontrivial and interesting proper- 
ties of the three-point function. Noting that we can always make A^ ' be zero by suitably 
redefining the field, the leading term that remains at late time is given by B U1 (A^) pro- 
vided neither _B m nor vanishes. Since this combination (1, 1, 1) provides the lowest 
decaying rate of the correlator, this is the only leading term that survives at late time. 
Therefore, Eq. (j34j) in this case becomes 

(0(t,£i)0(t,£ 2 )0(t,$O) = M 2 B in (A^) 3 (HR 12 y^(HR 23 y^(HR 31 y^ . (43) 
In this case, two-point function at late time becomes 

m,x x )^{t,x 2 )) =N\A^)\HR l2 y^. (44) 

Scaling behavior of Eqs. (I43p and (144)) coincide with the one given in [8HTT]. In these 
references, the scaling behavior was derived by combining the de Sitter isometries like we 
have done in this paper and the assumption that <fi (not at the level of the correlators) 
at sufficiently late time scales as (p ~ (—v) A an d this scaling directly enters the scaling 
of the correlator, for instance, {${j)i)(j>{j)2)) ~ ( — Vi) A (~ V2) A - On the other hand, as we 
have shown, the stochastic formalism can provide a necessary and sufficient condition in 
order for the above naive assumption to hold, which is given by A^B 1U ^ 0. In addition 
to this, the stochastic formalism also gives amplitudes of the correlators in terms of the 
calculable quantities. 
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What happens if Bni accidentally vanishes, which is possible for some models? In this 



case, the leading contribution to the three-point function comes from a term Bn 2 (A^) 2 A^ 



unless it vanishes fH Then, the two-point function is still given by Eq. (I44p . but the three- 



point function becomes 



(0(t, x^t, x 2 )<p(t, x 3 )> = Af 2 B 112 (AU) 2 AW 

(45) 

This is very different from Eq. ( 143|) on two counts; it does not obey the single power-law, 
and it is not given by a single term but by three terms that are mutually related by 
permutations. As far as we know, this type of three-point functions has been overlooked 
in literature. If B n2 (A^) A^ vanishes too, then we need to consider B 122 A^ (A^) or 
Bn-jA^ (A^) , whichever yields the lower decaying rate. If both of these have the same 
decaying rate, we must keep both of them in the correlator. Obviously, any of these leads 
to the multi-scaling expressions of the three-point function like Eq. (|45|) . These examples 
suggest that knowledge of the scaling behavior of the two-point function is not necessarily 
enough to know the scaling behavior of the three-point function. In principle, three-point 
function can exhibit more complicated structure than the native expectation. 

The single field with Bm ^ gives the following form of the four-point function, 

\n 13 K 2i -K23-K14/ ~~ 

This expression is the same as the one given in [10]. The function f uu for some limiting 
cases are given in Eqs. (138]) - fid]) . For instance, we find 

f un {z, z) = B 2 1U (A^)\-^, for z < 1. (47) 

Generally speaking, Bm = does not imply the vanishing of / mi but only changes the 
form of f l . Assuming that B\\ 2 does not vanish, / 1 corresponding to Eq. (147j) now 
becomes 

f lin {z, z) = Bf 12 (A^)\^-^, for z < 1, (48) 

which contains information of the second excited state $2- As it should be, the power of 
z is higher by (A 2 -AJ/H than Eq. (H7j) . 
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2A±— A2 A 2 A 2 

(HR 12 ) —(HR 23 y^(HR 31 y^ +2 perms. 




# 1 If the potential has a reflection symmetry V(<p) = V (—</>), &i(4>) for the bound state is an odd 
function, so that Bm in Eq. pip vanishes. In this case, however, -B112 should also vanish because of the 
same symmetry. 
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